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Abstract 

B c meson production in Pb+Pb collisions at y/s = 2.76 A TcV is surveyed 
in both a statistical coalescence model and a transport model. The nuclear 
modification factor Raa is predicted to be between 2 and 18 in the most central 
collisions, which can help to confirm the regeneration mechanism. Momentum 
dependence is also investigated in the transport model. A strong suppression 
of the transverse momentum is found in central collisions accompanying the 
enhancement in yield. The spectrum and elliptic flow of B c are also discussed. 



1. Introduction 

Several kinds of signals of the hot quark matter created in relativistic heavy 
ion collisions are widely discussed in the literature. One of them is based on 
the hard probes, for example quenching jl| of the jets which pass through the 
hot medium. Another one is focused on the low energy particles, for example 
the collective flow of low px particles [H, which carries the information on 
hydrodynamical properties of the hot matter at the initial stage of the fireball 
evolution. 

One more significant idea utilizes heavy quarkonia as hard probes of the 
deconfinemcnt. The anomalous J/ip suppression in quark-gluon plasma was 
predicted theoretically [H a quarter of century ago. Such a suppression was 
indeed observed in heavy ion collisions Q. It was however suggested that this 
suppression can bc as well explained within the models that do not assume 
the quark deconfinement Moreover, it was suggested that not only dis- 

sociation of quarkonia but also the opposite process, recombination of heavy 
quark-antiquark pairs, can take place in the deconfined medium (lp| - [l9j . If the 
regeneration of quarkonia indeed takes place, it becomes especially significant 
with increasing collision energy since the number of heavy quarks per collision 
becomes larger. When the regeneration is dominate, quarkonia can be used 
to detect the early fireball through its spectrum and flow like the soft probes. 
Even enhancement of quarkonia production has been expected [l4j , which could 
be regarded as a direct evidence of regeneration mechanism. However, such an 
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enhancement of heavy quarkonia is not supported at RHIC energy by the sta- 



tistical model [ll|,l2Cj and has never been observed at RHIC or LHC [2l|-|26|. 

How can we verify the regeneration mechanism firmly? 

B c meson was found by CDF in 1998 [27| . Similar to heavy quarkonia it 
consists of a heavy quark and a heavy antiquark. But unlike in J/ip or T 
mesons, the quark and antiquark have different flavors. This has a drastic effect 
on the production cross section of the B c meson in elementary hadron-hadron 
collisions. Indeed, creation of one cc or bb in an elementary collision is sufficient 
for production of, respectively, J/ip or T. In contrast, production of a B c meson 
in hadron-hadron collisions requires creation of at least two heavy quark pairs, 
cc and bb, in the same collision. 

Unlike in elementary hadron-hadron reactions, the B c meson production in 
nucleus- nucleus collisions does not require production of two heavy quark pairs 
in the same elementary collision. Heavy quarks and antiquarks originated from 
different nucleon-nucleon collisions can recombine and form a B c meson. If such 
regeneration is possible, it leads significant enhancement of B c production in 
nucleus-nucleus collisions relative to proton-proton ones. 

It was suggested that, due to the regeneration, the B c meson can be ob- 
served in Au+Au collisions at RHIC, although its observation in proton-proton 
collisions at the same energy is hardly possible [28|. This prediction is still to 
be verified in experiments. 

The high beam energy of LHC make it possible to measure the B c mesons 
both in p+p collisions and in Pb+Pb collisions. The cross section of bottom in 
proton-proton collisions is about one order of magnitude smaller than that of 
charm, while the cross sections of B c and T are about 2 orders less than J/ip. 
For a very rough estimation, the yield of a meson in the regeneration process 
can be assumed to be proportional to the yield of each of its constituent quarks. 
Therefore if the regeneration is present in nucleus-nucleus collisions, the yield of 
T mesons would remains two order of magnitude smaller than that of J/ip, while 
the multiplicity of B c mesons should be only one order of magnitude smaller in 
comparison to J/ip. Thus the nuclear modification factor Raa f° r B c can be 
roughly one order of magnitude larger than that of J/ip and T, that is Raa > 1 
which implies enhancement of B c instead of suppression. If such an effect is 
found, it will be a firm confirmation of the regeneration. 

In this paper, we will calculate the yield of B c meson in central Pb+Pb 
collisions in the statistical coalescence model, and then discuss more properties 
including the momentum dependence in a detailed transport model. We take 
7i = c = £;B = Iin the following. 



2. Baseline 

The B c meson is similar to quarkonia in the sense that they are composed 
of heavy quarks, so that the non-relativistic approximation can bc applied to 
describe the interaction between them. On the other hand, the decay modes 
of B c is totally different from quarkonia. For the excited B c states with the 
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mass below the B + D threshold, the only decay mode, except weak decay, 
is feeding down to the lower states until the ground state is reached, due to 
the conservation of charm and bottom. Thus the contribution from excited 
states are even more significant than in the case of quarkonia. The properties of 
ground and excited B c mesons are discussed in different potential models (29j . 
In order to be simple and to generalize it to finite temperature later in the 
transport model, we calculate the spectrum of B c by solving the Schrodingcr 
equation with a Cornell potential V(r) = —a/r + err. By this method, the 
mass spectra of charmonia and bottomonia can be well reproduced with the 
parameters a = tt/12, a = 0.2 GeV 2 , m c = 1.25 GeV, and m b = 4.25 GeV do|. 
With exactly the same parameters, we can obtain the mass spectrum of B c below 
the threshold as m Bc (lS) = 6.36 GeV, m Be {lP) = 6.72 GeV, m Bc (2S) = 6.90 
GeV, and m,B a (lD) — 6.98 GeV. Since the binding energy of the ID state in 
vacuum is only about 170 MeV^ T c , we neglect the contribution of this state 
in the following for simplicity. The contribution of this state in the statistical 
coalescence model is only a few percent. Note that there is already experimental 
data of B+ with quantum number J p — 0~ and mass m = 6.28 GeV [3l[, we 
use the experimental value instead in the following for the 15 state. 

The cross section of B+ at y/s = 1.96 TeV is measured by CDF as da/dy ( pr > 
6 GeV) = 15.5 ± 5.0 nb with rapidity \y\ < 1 According to PYTHIA 
the cross section for pr > at LHC energy y/s = 2.76 TeV is larger by 
the factor of about 4. Thus we obtain the inclusive cross section of B c as 
da/dy(pr > 0) = 62 nb. The cross section of charm and bottom are estimated 



as dcr c5 /dy\ y=Q = 620 /ib [3J, |35j and da b (,/dy\ y=0 = 20 (jh |36|, respectively. 
The elliptic flow and high px suppression of open heavy flavors implies strong 
interaction between heavy quarks and the fireball [37l ] , therefore we simply take 
the thermal momentum distribution for the heavy quarks. 

In the transport model, the momentum distribution of the initially produced 
B c s is required. It is parameterized in the power-law form 

2(n-l) _ , p\ Y" d<T (!) 



dydpj 



V J PT) {n-2)(pl) PT { 1 ' (n-2)G£) 



y=0 



with n = 4.16 and = 25.1 GeV 2 estimated from PYTHIA [H] simulation. 
To consider the reaction rates of different B c states, the cross section and the 
branch ratios of each of them are necessary. Since all the excited states feed down 
to the ground state without weak decays considered, we take the branch ratio as 
100%. The ratio of direct production cross section of ip' to J/ip is around 0.3 at 
low pt dil . As an analog, we take the cross section of each excited component 
of B c as 0.3 times of the ground state. That is a dlr (lS) : a dlr (lP) : a dir (2S) = 
10 : 9 : 3, where we have counted the degeneracy. To include the Cronin 
effect, a Gauss smearing is used to modify the initial momentum dependence 
so that the mean p T of the J/tjj in A+A collisions is larger than that in p+p 
collisions, that is (p T ) aa = (Pt)pp + a gN ' h where I is the total length of the 
path that the gluons pass through the nuclei before merging into a B c , with 
the broadening factor a g N taken as 0.2 GeV 2 /fm. Since the initially produced 
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B c s suffer strong suppression, the dependence of the final observations on the 
parameters discussed in this paragraph for the initial production is very weak. 



3. Statistical Coalescence Model 

The statistical hadronization model (SHM) and statistical coalescence model 
(SCM) have made grea t success in light and heavy hadrons respectively with 
few parameters [ill. |39|. This suggested the idea that hadronization of heavy 
hadrons can also be described within the statistical approach. In contrast to 
light quarks, the masses of c and b quarks are much larger than the typical 
temperature of the fireball. Therefore, production of cc and bb pairs at the 
thermal stage of the reaction can be neglected, though pre-equalibration can 
slightly alter the picture for cc production at LHC energy [40[. Practically 
all heavy flavor pairs are produced at the initial stage of the nucleus-nucleus 
reaction in hard parton collisions. It is assumed that their numbers remains 
approximately unchanged during the fireball evolution. 

The basic idea of SCM for heavy flavors is expressed in the following balance 



equation [11 



N QQ = \9Q N oQ+9 2 Q N% Q , (2) 

where Nqq is the number of directly produced heavy quarks, N*q and N^q are 
numbers of open and hidden heavy flavors for the hadron gas in complete (in- 
cluding heavy flavors) chemical equilibrium and zero value of the corresponding 
heavy flavor chemical potential. A fugacity gQ is introduced to describe both 
the conservation of Q and that of Q. When Nqq is small, the event by event 
fluctuation becomes important. This effect can be included by a modification 
factor (1 + 1 /Nqq) in the hidden part, that is N q q = (1 + 1/ N qq)9q N qq, as 
long as the heavy flavor is mainly open heavy hadrons (l2| . This modification 



leads to similar results of the widely used canonical ensemble modification [41 
with a deviation within several percent. In order to take the B c meson into 
consideration, we generalize the above equation into the following form 

N% r = \g c XK h c +g b N%)+[l + ^jglNlt (3) 

Nt = ^6«+ 9c <J+||l + ^U 6 Xt (4) 

where N oc and N a b do not include B c . This formula can easily be generalized 
for hadrons with even more heavy quarks. The leading order approximation for 
B c is simply 



... m^Njn 

Nb ° Nthmn — ( 5 ) 

oh oc 
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Higher orders counts for the strict conservation of heavy quarks and the balance 
between open and hidden heavy flavors. In the following calculation, the full 
form (J3J31) is solved. 

The temperature and volume of the fireball are taken as T = 164 McV 
and Va v =i = 4160 fm 3 respectively [42| for central Pb+Pb collisions at y/s = 
2.76 A TeV, and N c5 , N b i are derived from Glauber model with the cross sections 
discussed previously. All the listed charmed and bottomed hadrons with a mass 
m and spin J in the particle list from Particle Data Group [3l| are included. 
The results are as follows. g c = 31.1, g b = 2.39 x 10 8 , N ob /(2N^ r ) = 98.3%, 
N B+/ N bi r = L08% ; N r/N^ r = 0.60%. Note that 



Nj c r of; _ 62 nb 

rdir —dir 



0.31%. (6) 
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OO 00 r 

The nuclear modification factor can be calculated as 

R N Bc N Bc /N^ _1.08%_ 

AA iV|«" N^/N^f 0.31% ' y ' 

Only the scalar component of B c mesons has been observed in experiments. If 
we take the probably existent vector state of B c and the isospin partner of B* 
into consideration, the results are modified as follows. g c = 31.1, gt = 1.83 x 10 8 , 
N ob /(2N^ r ) = 96.4%, N B +/N^ r = 3.30%, N Ts /N^ r = 0.35%, which results 
in Raa = 3.30/0.31 = 11. Further more, if B C (1P) and B C (2S) states are also 
considered, the Raa can be as large as 13. Thus from the SCM, we do expect 
an enhancement of B c mesons in the central Pb+Pb collision at LHC energy 
with the Raa = 3.5 ~ 13. 



4. Transport Model 

In contrast to the statistical coalescence model, that assumes statistical 
equilibration of the heavy quark distribution among the hadrons at chemical 
freeze-out, the transport model takes into account deviation from the statistical 
equilibrium and trace the whole evolution of the fireball. In this model, the 
distribution of B c in the phase space is described by the function /s e (x, p, t) 
satisfying the transport equation 

(fit+v V)/ Bo =-af Bc +/3, (8) 

where v = p/Es c is the velocity of B c . The dissociation rate can be expressed 
as 

a( X ,p,t) = — [ , ^ k^fl h a d j F ~ T i , (9) 
1 ' P ' ' E B J {2irfE k ^ Jg 6{T d - T) K ' 

where fg h (k, u, T) is the thermal distribution of gluons, a(k,p) is the cross sec- 
tion of the gluon dissociation process + g b + c, which is obtained by 
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replacing the heavy quark mass ttiq with twice the reduced mass 2[i and substi- 
tuting the proper binding energy in the cross section for quarkonia obtained by 



OPE method [281. |43H45| . In B c rest frame, it can be expressed as [281. 143H4 

(w/eis - 1) 3/2 



0isM = ^o- 



<7 1P (U>) = AA 



02sM = 16A 



(w/eis) 5 

(lj/6 1P - l) 1 ' 2 ^/^) 2 - 20(^/e 1P ) + 12) 
(wAip) 7 
( W /e 2S -l) 3 / 2 (w/e 2 s-3) 2 



(uj/e 2 s) 7 



where u is the gluon energy, A = 2 11 7r3 _3 (2^) _3 / 2 e 1 ^ 2 , \i = (m&m c ) / '(mj + 
m c ) is the reduced mass and ej is the binding energy of the B c state J. Here we 
use vacuum values, that is the binding energy is found as a difference between 
the sum masses of ground state B and D meson and the mass of corresponding 
B c state. To take into account the recoil effect due to the finite mass of B c , we 
further replace the binding energy e by the threshold ujq = (1 + e/(2m,B e )) e as 
in 0. 

The lifetime of B c with account for the gluon dissociation is shown in Fig. [TJ 
The binding energy of the ground state is much larger than that of the excited 
states, therefore the ground state lives longer. Comparing the lifetime of B c me- 
son to the radius of the nucleus R(Pb) ~ 6 fin, one finds that the ground state 
suffer strong suppression at the temperature above 400 McV and little suppres- 
sion below 300 McV, while the suppression to the excited states is already very 
strong at 200 MeV. 

The thcta functions in the numerator and denominator in ((9]) restrict the 
above process between the critical temperature T c of light hadrons and the 
dissociation temperature Td of B c mesons. Below T c , there is no gluons in the 
fireball. Above Td, the bottom and charm quarks are strongly screened, and 
can no longer form a bound state. 

With the heavy quark potential at finite temperature obtained from lattice 
QCD [13], the dissociation temperature Td of B c mesons can be calculated by 
solving the Schrodinger equation as in [30(. The results are listed in Table [TJ 
The dissociation temperature is calculated for two extreme cases. In the first 
one, the potential is assumed to be equal to the internal energy: V = U, and 
in the other one it is equal to the free energy V = F. Note that the difference 
between them is the entropy term. The former corresponds to a quick adiabatic 
screening, while the latter assumes a strong heat exchange between the heavy 
quark system and the medium (48j . Since the difference is remarkable, we will 
do calculations in both limits in parallel. 

Besides the dissociation of B c , the inverse process b + c — > B c + g is also 
considered in ([8]), which lies in the regeneration rate 

" = lk I (2^(2^p^' l ''' l ''' i(1+ ' ,)(2 *'' ,l ' ,|, ' + ' _ '' _s1, 
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Figure 1: Lifetime of B c meson with account for the gluon dissociation in QGP. The velocity 
of the QGP is taken as v = in this figure. 



States of B c 


IS 


IP 


25 


T d /T c (V = U) 


3.27 


1.59 


1.41 


T d /T c (V = F) 


1.51 







Table 1: The dissociation temperature of B c mesons scaled by the critical temperature T c . 
The heavy quark potential V is considered as the internal energy U and the free energy F at 
finite temperature, respectively. The dash means the meson does not survive above T c . 

where p, k, q c , are momenta of -B c , gluon, charm and bottom respectively. 
W(s) is the transition probability, which is related to that of the dissociation 
process by detailed balance of the cross section entering in Eq.©. f c , fi and 
f g are distribution functions of c and b quarks and gluons. f g is taken as 
thermal distribution as in the dissociation rate a. For simplicity, we also assume 
kinetic thermalization for the heavy quarks, and neglect the Pauli blocking 
effect. That is the distribution function /q = PqSq, where the momentum 
part /q is normalized Boltzmann distribution, and the density of heavy quark 
Pq is determined by the conservation law of heavy quarks 

dt P Q + V-(pQv) = 0, (10) 

where v is taken as the velocity of the medium. The initial condition is deter- 
mined by Glauber model. The regeneration process takes place in the temper- 
ature window T c < T < Tj. When T <T Cl There are no partons, therefore the 
previous process does not exists. When T > T^, the dissociation rate a is infin- 
ity, and the B c gets dissociated as soon as it forms, therefore the regeneration 
does not work in high temperature, either. Since the dissociation temperature 
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Td depends strongly on the form of the heavy quark potential, one would expect 
that the yield of B c from the regeneration process also depends strongly on the 
potential. 

As a background of the B c suppression and regeneration, the fireball is de- 
scribed by 2+1 dimensional ideal hydrodynamics with the assumption of boost 
invariance [49[. The equations of state (EoS) are taken as massive ideal gas of 
hadrons and partons in the confined and deconfined phases respectively with 
a first order phase transition at T c = 165 MeV as it was done for J/ip and 
T [Hoj . Multiplying this critical temperature by the ratio Td/T c in Table [TJ the 
dissociation temperature Td is obtained, regardless of the absolute value of T c in 
lattice simulation. Otherwise there could be even more enhancement of B c . The 
initial condition is decided by Glauber model, and the maximum temperature 
of the fireball in the most central collision is 485 MeV at the thermalization 
time To = 0.6 fm. 

The ccntrality dependence of the nuclear modification factor Raa in both 
limits of V = U and V = F are shown in Fig. [5] If the regeneration is ignored, 
the Raa of the initially produced B c is 0.25 for V = U and 0.01 for V = F, which 
are reasonable when compared to the experimental results of T. The Raa& of 
T(15 +2S + 3S) at RHIC |j| and T(15) at LHC (EH in central collisions are 
consistent with the simple assumption that the ground state survives while all 
the excited states melt in the fireball. One can expect that the suppression of 
B c is stronger than that of T. In the V = U limit, the maximum temperature 
is still lower than Td of B C (1S), and about half of the B C (IS) survives the 
gluon dissociation, which is less than that of T(15). In the V = F limit, 
the dissociation temperature is much lower, and the ground state suffers even 
stronger suppression. 

In both cases, the contribution from the initial production is two orders of 
magnitude smaller than that from the regeneration, and thus the population 
of B c is dominated by the latter. Note that a larger Raa of B c than that of 
T already implies the regeneration mechanism in the fireball. In both limits 
of our calculations, the Raa = 2.3 ~ 17.5 is over unity. The temperature for 
regeneration in the transport model is higher than that in SCM, thus one can 
not expect whether the yield of B c is smaller in the transport model or in the 
SCM. Actually, this interval covers the SCM results and also implies strong 
enhancement. The remarkable difference between the two limits comes from 
two factors. First, the excited states does not survive at the free energy limit, 
while they play an important role in the internal energy limit. In the most 
central collisions, the contribution from the ground state and the exited states 
are almost the same for V = U. Second, the dissociation temperature of the 
ground state is much higher in the internal energy limit, which allows more 
regeneration in the hot fireball. 

When the regeneration inside the fireball is confirmed, the nuclear modifi- 
cation factor loses its meaning as a survival probability, because the observed 
B c mesons are mostly regenerated instead of surviving. Thus the Raa depends 
not only on how the B c interacts with the medium, but also on the production 
cross section of B c and heavy quarks in hard nucleoli- nucleon collisions. 




Figure 2: The nuclear modification factor Raa of B c meson as a function of the participant 
nucleon number Afp ar ^ in V = U (left) and V = -F(right) limits. The long-dashed, dashed, 
and solid curves are the contribution from initial production, regeneration and the total re- 
spectively. 



To characterize the producing and/or surviving ability of B c from the regen- 
eration, we take a ratio of the final yield of B c mesons over those of the heavy 
quarks in unit rapidity region [52j as 

dN B Jdy 

dN c /dydN b /dy' K ' 

While Raa compares the final yield of B c to the initial yield, F compares 
it to the source of regeneration, and thus becomes independent of the initial 
cross sections of B c and heavy quarks when the regeneration dominates. The 
final ratio F at LHC energy is shown in Fig. [3] compared with that of the 
regenerated B c at RHIC. The final ratio is larger in the V = U case, as one 
would expect from the dissociation temperatures. Since the temperature at 
LHC is much higher, the final ratio is smaller than the corresponding case at 
RHIC, which is consistent with the picture of color screening. The smaller F in 
central collisions also implies smaller producing and/or surviving ability of B c 
in hot large medium, and the growth of Raa with N-pg^ is mainly due to the 
increasing multiplicity of heavy quarks. Meanwhile the yield is very sensitive to 
the potential as we found in Raa, so that the F with V = F at LHC is even 
smaller than that with V = U at RHIC. The different behaviors between RHIC 
and LHC at small -?Vp ar t is artificial. That is because the initial production 
is not included in calculations for RHIC in Fig. G2 and F as a measurement 
of producing and/or surving ability only makes sense when the regeneration 
dominates. 

The centrality dependence of F in SCM is similar, as shown in FigJU In 
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Figure 3: The final ratio ^(defined in 11H 1 of B c meson as a function of the participant nucleoli 
number A r p ar ^. The long-dashed and solid lines are V = F and V = U limit respectively. 
The black thick lines are for LHC energy, while the red thin lines are for the regenerated B c 
mesons at RHIC energy as a comparison. (Color on line) 



order to apply the SCM model in different centrality, we assumed that the 
volume of the fireball is proportional to the multiplicity of charged particles, 



which is measured by experiments |53|- The yield of heavy quarks is from a 



Glauber model as before. The result can be understood from ([5]), which implies 

4N th N 1 ! 1 An th n tl l 1 1 

p _ ^ iv oc iv ob _ ^' L oc n ob 1 J_ 

Nf e nf c V V 

The volume increases with iVpaxti thus the final ratio F decreases, which implies 
the fact that when the system is more dilute, it is more difficult for the heavy 
quarks to meet and combine with each other. When considering F ■ V, there is 
almost N p scaling, with a deviation within a few percent in peripheral collisions. 
The deviation comes from the change in heavy quark density and the event by 
event fluctuation. 

One of the advantages of the transport model is to investigate the momentum 
of particles, which is sensitive to the dynamics. The initially produced B c 
mesons come from the hard collisions, while the regenerated B c s are merged 
from the heavy quarks that are softened by the medium. The typical energy 
scale of the medium is the temperature, which is much smaller than the typical 
energies of hard protons in the initial beams. Thus the transverse momentum 
carried by the initial production is obviously larger than that of the regeneration. 
When the regeneration becomes important, there is a suppression in (p%) , as 
shown in Fig.O The same phenomenon is also expected for J/ip at LHC [54| . 



10 



0.04 



0.03 
F 0.02 



0.01 



0.00 

100 200 300 400 
Npart 



Figure 4: The final ratio ^(defined in {TTJ) of B c meson as a function of the participant 
nucleon number Np ar ^ from SCM calculation. The red long-dashed and black solid lines are 
calculations of observed scalar B c state only and that including excited B c states (as listed 
in Table [TJ, respectively. (Color on line) 
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Figure 5: The average transverse momentum square (pj,) as a function of the participant 
nucleon number Afp ar ^ in V = U (left) and V = -F(right) limits. The long-dashed, dashed, 
and solid curves are for initial production, regeneration and the total, respectively. 

Since thermalized heavy particles from a point-like thermal source follows 
Boltzmann distribution, that is dN/dp oc exp(— Et/T) at mid-rapidity, we plot 
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\n(dN / dp) at mid-rapidity as a function of Et in Fig. [B] for central Pb+Pb col- 
lisions. It would be a straight line for Boltzmann distribution. The softening of 
the spectrum in Pb+Pb collisions is obvious through the change of the slope. If 
a superposition of Boltzmann distributions at different temperatures is consid- 
ered, only a concave curve is expected. However, in the internal energy limit, 
it is a convex curve at Et ~ 9 GcV. Such a behavior mainly attributes to the 
suppression of excited states. If we switch off the gluon dissociation process of 
the excited states, the obvious bend at low pt disappears. The effective tem- 
perature extracted from the spectrum at low pr according to the Boltzmann 
distribution is about 550 McV, which is above the maximum temperature input 
of the fireball. This mainly results from the radial flow. When the velocity of 
the fireball is switched off in the regeneration process, the effective temperatures 
are around 250 MeV, while a calculation of average temperature directly from 
the transport model is about 230 MeV. Since the radial flow seems large, one 
would also expect large elliptic flow. 
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Figure 6: B c meson spectrum with respect to the transverse energy Et = \J Vr + m2 ■ I n 
V = [/(left) and V = F(right) limits. The long-dashed, dashed, and solid curves are for initial 
production, regeneration and the total, respectively. 

The elliptic flow at b = 8.4 fm is shown in Fig. [7] The initially produced 
B c s do not thermalizcd with the medium, and thus carries small elliptic flow. 
The non-zero flow comes from the suppression process. The regenerated B c s 
are born inside the fireball, thus inherit a relatively large flow of the medium 
through heavy quarks. The total is dominated by the regeneration at low px, 
and decrease at high pt when the initial production becomes important. Since 
ideal hydrodynamics and the kinetic thermalization of the heavy quarks are 
assumed, the flow of regenerated B c at high pt is not reliable. Note that the 
assumption of the thermalization of the momenta of bottom quarks is essential 
in this i>2 calculation. If the momenta of bottom quarks is not thermalized, the 
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flow of B c could be smaller, which is to be verified in the experiments in the 
future. 
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Figure 7: The elliptic flow V2 of B c mesons as a function of the transverse momentum px at 
impact parameter b = 8.4 fm in Pb+Pb collisions. Two limits of the heavy quark potential 
V = [/(left) and V = .F(right) are shown. The long-dashed, dashed, and solid curves are for 
initial production, regeneration and the total, respectively. 



5. Conclusion 

Based on the results of both the statistical coalescence model and the trans- 
port model, an enhancement of the B c mesons yield in Pb+Pb collisions rela- 
tive to p+p collisions at yfs = 2.76 TcV is predicted. If such an enhancement 
is observed, we can conclude firmly that the regeneration of B c in the fireball 
happens, thus the J/ip production in Pb+Pb collisions at LHC is most likely 
dominated by the regeneration mechanism, and new variables such as F in ((4]) 
can be introduced to characterize the properties of quarkonia in the fireball. In 
the opposite case, i.e. if no B c enhancement is observed, one has to conclude 
that the observed T yield is dominated by the initial production. According to 
our transport results, the transverse momentum of B c is suppressed accompa- 
nying the enhancement in yield, since the regenerated B c s are produced at a 
lower energy scale than the initial production. The pr spectrum is much softer 
than that in pp collisions, and in the V = U limit, there is a suppression at low 
Pt compared to the thermal distribution, which results from the suppression of 
the excited states. Hopefully, the study of B c mesons will shed more light on 
the production and suppression of heavy quarkonia and properties of the hot 
quark matter. 
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